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Abstract 

We solve a two dimensional model for polymer chain folding in the presence of mechanical pulling 
force (/) exactly using equilibrium statistical mechanics. Using analytically derived expression for the 
partition function we determine the phase diagram for the model in the /-temperature (T) plane. A 
square root singularity in the susceptibility indicates a second order phase transition from a folded to 
an unfolded state at a critical force (/ c ) in the thermodynamic limit of infinitely long polymer chain. 
Surprisingly, the temperature dependence of f c shows a reentrant phase transition, which is reflected 
in an increase in f c as T increases below a threshold value. For a range of / values, the unfolded 
state is stable at both low and high temperatures. The high temperature unfolded state is stabilized 
by entropy whereas the low temperature unfolded state is dominated by favorable energy. The exact 
calculation could serve as a bench mark for testing approximate theories that are used in analyzing 
single molecule pulling experiments. 



Since the pioneering experiments [TJ? , [2] demonstrating that titin, a large protein with mul- 
tiple subunits, can be unfolded upon application of mechanical force, /, single molecule pulling 
experiments with increasing sophistication have been used to extract the folding landscapes of 
proteins [U [5] and RNA [61 [7] . The experimentally accessible coordinate in these experiments 
is the molecular extension, R, which is conjugate to / and represents typically the distance be- 
tween the ends of the molecule of interest. Parameters describing the folding landscape, such as 
the barrier to unfolding at / = 0, stability of the bound or folded states with respect to the un- 
bound states, and location of the transition state are extracted from measurements (unbinding 
rates as a function of / for example) by assuming that a single reaction coordinate, R, suffices. 
Although such measurements and interpretations have provided insights into some aspects of 
the folding mechanisms the inability to confirm the accuracy of the extracted parameters by us- 
ing independent measurements of the key quantities remains a major, but not widely discussed, 
problem. This situation is exacerbated by a paucity of exactly solvable models (however see 
[8]), which can be used to test the accuracy of various commonly used approximate procedures 
to extract the folding landscape parameters from measured trajectories. 

Motivated in part for the reasons stated above, in this work we obtain an exact solution 
on the effect of force on a remarkable toy model introduced over forty years ago by Zwanzig 
and Lauritzen (ZL) [H] in the context of polymer crystallization. We construct the equilibrium 
phases of the ZL model in the [f,T] variables where T is the temperature. Previously such phase 
diagrams have been obtained numerically using models introduced to understand protein folding 
\10\ ITT] and also by theoretical methods for lattice homopolymer models [12]. As stated above, 
we were motivated to undertake this work in part to discover models for which exact solutions 
can be found so that the accuracy of approximate methods can be assessed. In addition, there 
is growing interest in the collapse transition of self attracting polymer chains, because of the 
relevance to protein folding [T3HT5] . Of course, chain folding is also at the heart of self-assembly 
of proteins and RNA. Such transitions are caused by a competition between the self attraction 
between segments of the polymer (tending to collapse the chain) and the conformational entropy 
that favors expansion. Typically, at low temperatures the collapsed state is stable dominated 
by energy or enthalpy. Increase in T populates extended conformations that are stabilized by 
entropy. Since, the ZL model includes both these effects we suspect that certain generic aspects 
of chain folding can also emerge from a deeper study of this interesting model. Although the ZL 
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model or the extension investigated here and elsewhere [TBI ITT] may not be directly applicable 
to proteins the study of such models is interesting in its own right, and for clarifying unusual 
aspects of phase transitions in polymers [18] . 

The ZL model was introduced as a caricature of polymer crystallization, which in the absence 
of / exhibits a phase transition that can be treated exactly by equilibrium statistical mechanics. 
The ZL model leads to a novel second order-like phase transition from an extended state to 
chain folded state in the thermodynamic limit of infinite chain length. The exact solution to 
the ZL model provided here shows folding-unfolding transition in the presence of mechanical 
force (/). In particular, we show that a second order phase transition also occurs from a folded 
state to an extended state at a critical force f c at a fixed T. The exact calculation of the phase 
boundary separating the folded and extended states show a reentrant phase transition in which 
the extended state is stable at both high and low temperatures at finite /. 

Zwanzig-Lauritzen model: 

The ZL model considers a long polymer molecule whose conformations are restricted to a 
plane. At / = the chain molecule folds into lamellar structures at low temperatures reminiscent 
of conformations adopted by crystalizing polymers. A conformation of the ZL chain, displayed 
in Fig. 1, shows that it can fold upon itself by paying a pending penalty, which is compensated 
by attraction between the chain segments. Consequently, at / = the chain molecule folds into 
lamellar structures at low temperatures. The molecule folds into N segments ranging from one 
to infinity, 



The folds are uniform occupying a length q (Fig. 1). Let the length of i th segment be X{. It can 
have any length, 



1 < N < oo. 



(1) 



< Xi < oo, 




with the constraint that the total length L, 



N 



J2xi + (N -l)q = L 



(3) 



i=l 



be fixed. 
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We apply a mechanical force at one end of the molecule, say /, with the other end fixed. The 
total energy of a conformation is associated with (1) a bending energy u for each fold (Fig. 1), 
(2) the interaction energy between any pair of neighboring segments of length X{ and Xj+i, 
and (3) stretching energy due mechanical force at one end. Two neighboring segments interact 
with short-ranged van der Waals forces. ZL assumed that the attractive energy of interaction 
between i t h and (i + l) t h segments is — J mm(x„Xj + i), where min(a, b) stands for the smaller 
of the length of the two segments. The proportionality constant J is taken to be positive. The 

— * 

stretching energy due to the applied force / is taken to be, 

- / • (r) (4) 

where, f is the end-to-end distance of the molecule. In the present calculations we assume that 
the force is applied along the x direction. The total energy of a given conformation of the chain 
molecule is 

N-l 

E(xi,....x n ) = (N — l)u — J min(xi, Xi+i) (5) 

i=l 

-f x (N-l)q-f y (x N +p) 

where (N — l)q is the end-to-end distance along the x direction and xn + P is the end to end 
distance along y direction with arbitrary p, as shown in Fig. 1. 

At first glance the ZL model resembles a one-dimensional model with short range interactions, 
which cannot exhibit a phase transition. However, ZL pointed out that this model is analogous 
to the familiar one dimensional Ising model with N possible states for each of the N sites. A 
two dimensional Ising model can be thought of as a one dimensional Ising system in which each 
site has 2 N sites. Thus, apparently the ZL model has enough of the two dimensional features to 
exhibit an interesting second order transition from extended to folded states even in the absence 
of force. 

Evaluation of the partition function in the presence of force: 

We take the continuum limit of the model. In this limit, the canonical partition function 
Z(L), which is the sum of the Boltzmann factors e~ ,3EN of all configurations of a polymer 
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molecule subject to the constraint of total length L, is 

°° POO POO 

Z{L) = / ' ' ' / dxi ■ ■ ■ dxN (6) 

nr_i JO JO 



N=l 



N 

-PE N 



5(J2xi + (N-l)q-L)e 

4=1 

with (3 = j^f, and Kb is Boltzmann's constant. As we sketch below the partition function 
Z(L) can be evaluated analytically following the strategy developed by ZL. We start with Eq.(|6]) 
and perform a Laplace transform of the partition function with respect to length L to eliminate 
the delta function on the right hand side of equation Eq.Q. Then we invert the transform to 
obtain Z(L). Using the following relation, 

min(a, b) — -(a + b) — -\a — b\, (7) 

we can write the potential energy in the following form 



E 



N 



(N -l)u- -J^2( Xl + x l+1 ) + 



2 

1 



1 N-1 

-J ^2 \ X i ~ X i+l\) ~ f*( N - !)9 - fy( x N +p). 



We write Eq.([8j)as, 

E N 



2 

1 



- 1 

(JV-l)u- J22x i + -J(x 1 + x N (9) 
1 

JV-l 

+ ^2\xi- x i+1 \) - f x (N - l)q - f y (x N +p). 



The sum Xi in equation (JoJ) can be replaced by L — (N — l)q, leading to 



E N = (N-l)(u + Jq)-JL + ^J( Xl +x N (10) 



N-1 



+ ^2\xi- x i+1 \) - f x (N - l)q - f y (x N +p). 
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By substituting Eq.(lO) into Eq.^6j) we obtain, 

oo 

Z(L) = ^[e~ Ku+Jq ~ Uq) ] N ' 1 e pjL 

N=l 

poo poo 1 v 

/ ••' / dxi ■ ■ ■ cIxn8(S~2 Xj + (N — l)q — L) 

Jo Jo i=1 

1 1 N ^ 

)[--/3J(x 1 + x N ) + - J ^2 



11) 



N 



exp[ 



Taking the Laplace transform on both side to eliminate the delta function on the right hand 



+Jq-f*q)]N-l 



side of equation (11), we obtain, 

poo 00 

/ dLe-^ + ^ L Z(L) = J2[^ {u 

J** N=l 

roc poo 

e-^ N ~^ / ... / dx 1 ---dx N e~^* 
Jo Jo 

1 1 N ~ l 

exp[--j3J(xi + X N ) + - J ^ \ X i ~ 



(12) 



,Pfv(xN+p) 



The integrals in Eq.(12) are evaluated using the following method. We define a partial 



generating function with the initial point fixed at the origin. Let us define a recursion relation 
of the partial generating function with fixed x^ = x. 



g(z,x) = te 



g(x, y) + e ~ <x 



(13) 



where t — e P( u+Jq f*i) anc [ j = I {3 J. The needed generating function is 



G(t,T]) = / e^ x g(z = 0,x). 
Jo 



(14) 



The term J °° g(x,y) in equation (13) is the generating function in the absence of mechanical 
force, which is explicitly calculated in [H] by solving the differential equation given by, 

r]a J v {ae~ vx / 2 ) 



dyg(x,y) 



2t J„_i(a) 



(15) 
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13 



14 



and 
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with, a = (^f ) 1//2 and v = 27/77 and u is replaced by u + Jq — f x q. Using Eqs 
we can evaluate the generating function in presence of mechanical force exactly. We obtain 

<*■•>- 

ya" +1 P F q [{l - f + v\, {1 + v, 2 - f + u}; -f\ 
J,_ 1 (a)r(z/ + 1) 

with a = and 77 > {3f y — 27. Hence we get, 

POO 

/ dLe- (r?+/3J)i Z(L) = e^ u + Jq - f * q) e m G{t,i). (17) 
The partition function Z(L) can be found by taking the inverse Laplace transform, 

1 pc+ioo 

Z(L) = eP {u+Jq - f * q) - — -/ drie vL+m+ P JL G(t,ri). (18) 

The contour of integration is a straight line, parallel to the imaginary axis, to the right of all 
singularities of the integrand. 

We now analyze the analytic properties of G(t, rf) in the complex rj plane. Originally, the 
function G{rj) is calculated for real 77 with the imposed condition 77 — (3f y + 27 > 0. Now, we 
wish to analytically continue G(rj) off the real 77 axis. The Bessel function has branch points at 
the origin, and the usual cut along the negative real axis, but the related function, 

A^ 1 (z) = T(u)( Z -)-^J^ 1 (z) (19) 

is an entire function of z and a meromorphic function of order v — 1. We start with real positive 
77 and hence real positive u, where the function G is well defined, and can be written as, 

T(v)T(l - - + u) 

G »^)° 1X^1)1X2-^1) (20) 

r/a 2 p F g [{l-| + 4;{l + ^2-|};^] 



(2)i-"A I/ _ 1 (a) 

For simplicity, we assume that the pulling force is along x direction (Fig. 1). The function 
pF g [{l — § + v}, {1 + v, 2 — |}, = f-\ reduces to A u . Now A v and A u _i are analytic functions of 
complex v everywhere except at singularities of the gamma function. At these points, which are 
negative integers, they have simple pole. But, in the ratio, the singularities cancel out exactly 
at these points. This implies the singularities of G in the complex plane are determined solely 
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by the zeros of the denominator, which are now poles. This accomplishes the desired analytic 
continuation of the ratio of the positive real axis into the entire complex plane. 
Let us define, 

2 8 

° = (Jjyji exp [~ 2 + Jq ~ ^ 

so that, 

a = av. (22) 



In order to evaluate the contour integral in Eq. [L8j we need to evaluate the zeros of the 
denominator namely the zeros of the Bessel function equation, 

J v -x{av) = 0. (23) 

Each of the zeros will be function of o as well as the external force. To each zero, there 
corresponds to a particular value of rj, which are all real. The sequence of all zeros of J„_i(cn/), 
corresponds to the sequence of all real rj's, which is given by, 

Vi > V2 > V3 > ■ ■ ■ ■ (24) 

indicating the largest one is r/ 1 . Let Ri denote the residue at each pole of G(r]) corresponds to 
each 77 j. We obtain, 

Z(L) = exp(/3( M + JL + Jq- f x q)) ^ R 4 exp(L^). (25) 

i 

Phase Transition: 

We are interested in the thermodynamic limit, corresponding to the length of the polymer 
molecule going to infinity. In the large L limit, the largest 771 will dominate the sum. Conse- 
quently, the partition function becomes, 

lim [L~ l log Z{L)) = (3J + m(<r, (3, f x ). (26) 

L— >oo 

We have neglected terms that vanish in the large L limit. Thus, the evaluation of Z(L) as well 



as the free energy reduces to the computation of a particular zero of the equation (23), that 
gives the largest r\\. 

The behavior of a for the largest order is known to be, 

a = z/-l + a 1 (z/-l) 1/3 (27) 



for large z/,and 

a = av. (28) 

Comparing these equations, we conclude that for a < 1, there is no solution for positive v. 
There is an asymptotic solution v = oo, and for a > 1, there is a unique solution < v < oo. 
We assume here that the fold energy u is positive, a is monotonically increasing function of both 
temperature and force. Thus, a is unity at the transition point determined by either 

2 exp[-^(u + Jq-f x q)] = l, (29) 



for fixed f x , or 

2 exp[-^( M + Jq - f c q)} = 1, (30) 

for arbitrary /3. Equation (30) represents the exact form of the critical force for arbitrary (5. 

In the limit f3 — > j3 c from above or f x — > f c from above at a constant 0, v approaches oo, so 
that, 

av = v + au 1/3 . (31) 

It follows that 

v = a 3 /> - I)- 3 / 2 , (32) 
if T ~ T c + or f x ~ fc - or equivalently, 

u - 1 = a' 3 / 2 (a - if 2 . (33) 

Thermodynamics 

We discuss the thermodynamic properties of the model in the long chain limit. Let A denote 
the Helmholtz free energy of the system. In the thermodynamic limit, the free energy per unit 
length is given by, 

lim(y) = -J-/rV (34) 
By defining the reduced temperature T" = 4j, we obtain the reduced free energy, 

A' = hm ( A) = _l _ TV (35) 

The reduced mean energy is given by, 
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The average end-to-end distance in the direction of f x is < (N — l)q > can be found by 
differentiating the partition function with respect to f x , 

<(iV-l)g>=/3- 1 ^|T'. (37) 
dfx 

In the thermodynamic limit, the average end-to-end distance per unit length, I, along the x 
direction is, 

1 f)A' 
I = lim (j<{N- l)q >) = (38) 

L^oo L Of 

We are interested in the vicinity of the transition point, which is in the neighborhood of a = 1. 
Suppose that, a = 1, when X" = 2^(7) when f x is fixed, Taylor series expansion of a is given by, 

°{T') = 1 + [^r + U ~ 1 ^ 1 - Jfyl (T' ~T C ) + ... (39) 
Similarly, for arbitrary T", Taylor series expansion of a around the critical force f c is given by, 

* = l + [2^](/.-/c) + "- (4°) 



First consider the case with = 0. The expansion in (39) shows that v~ x is proportional to 
three halves power of the temperature deviation from T c . Now, for T' < T c , a is less than unity 
and 771 = 0. When T' becomes T c + , a becomes greater than unity. The free energy is given by, 



if T' > T c , and, 



L 2T C 2T C 2 



(41) 



if T < T c . 

Similarly, for specific heat, 

if T' < T r , and 



A' = -1 (42) 



, 9E . 

^)t = (43) 



= \[T C \[^ + ^^] 3/2 (T' - T c )-W (44) 

if T" > T c . Specific heat shows inverse square singularity at the transition temperature. The 
discontinuity in specific heat shows a second order phase transition at T' = T c for fixed f x =0. 

It can be easily seen that in the absence of force, and for negative value of the bending 
energy u, in the interval — 4- — Jq < u < —Jq, a is infinite for both T' — and T' = 00, with 
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a single minimum below unity at a finite temperature. It follows that a will intersect unity at 
two temperatures and hence there are two phase transitions. At low enough temperatures, with 
negative bending energy, increasing the number of folds will minimize the total energy. This 
leads to the stabilization of the stable folded state. Such a possibility is not realized for positive 
fold energy. 



Now, we fix T and vary the applied mechanical force. The expansion in (40) shows that v~ x 
is proportional to the three halves power of the force deviation from f c . For f x < f c , a is less 
than unity and r\\ = 0. As f x — > /+ , a becomes greater than unity. The free energy is given by 

A' = -1 (45) 

if f x < fc and 

A' = -l-[^77] 3/ U-/ c ) 3/2 (46) 

if fx > U 

The average end to end distance I per unit length in the thermodynamic limit is found in 



(38). If a < 1, then we have 771 = and the mean distance per unit length is zero. In the 
thermodynamic limit, any finite end-to-end distance corresponds to I = 0. The average end- 
to-end distance per unit length vanishes when f x < f c corresponding to the folded state, and 
increases with the square root of the force difference above f c for arbitrary T' representing the 
unfolded state as follows, 

/ = f x <f c , (47) 



and 



I ] 3/2 (/*-/c) 1/2 , (4* 

2 [ 2a 1 T'J i KJ J ' ' V 



if, f x > f c . Thus, the folded-unfolded transition in this model is second order even in the 

dl 



presence of force. Our calculations further reveal that the susceptibility \ = wr °f the polymer 



diverges as 

*«(/,- fc)- 1/2 (49) 

as f x approaches the transition point from above. The divergence in susceptibility at transition 
point shows the second order phase transition with the characteristic critical exponent 1/2. 
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Phase diagram in the [/, T] plane: 



In Fig. 2, we show the regions of / — T plane, where the chain is extended (white) and where 
it is chain folded (shaded) for this model. Remarkably, the phase diagram exhibits a reentrant 
behavior. The boundary separating the two phases occurs when 

2T /l/2 e -|^ e JM 7 = L ( 5Q ) 

This is an exact form for the temperature dependence of the critical force, f c . The transition 
curve separates the region of extended chain and chain folded configurations. Upon crossing the 
phase boundary from below the folded chain undergoes a transition to the extended state via 



second order transition. From Eq.(50), we can write 



f u 

J jq=( 1 + q)-T' log(4) - T' log(T'). (51) 
The transition line in Fig. 2 increases at low temperatures due to the leading term T' log T" 



in Eq.(51). In the interval u + Jq < f c q < u + Jq + j-, the folding transition occurs at a 



critical force over a range of temperatures. However, upon further lowering the temperature 
the chain becomes extended a process that is reminiscent of 'cold denaturation'. Hence, the 
extended state is stable at both high and low temperatures, implying that there is reentrant 
phase transition. At low enough temperatures, for the case of f c q > u + Jg, large end-to-end 
distance minimizes the total energy corresponding to a stable extended state. In particular, 
Fig. 2 shows the unusual behavior that at low enough temperatures the critical force required 
to unfold the chain is less that at a higher temperature. For example, f c at T = 0.05 is less than 
f c at T = 0.1. 



Conclusions: 



The exact solution of the ZL model for polymer chain folding in the presence of force shows a 
second order ordering transition from a folded to an unfolded state provided the values of force 
are less than a critical value, f c {T). This prediction is not in accord with mean field theory, which 
suggests that the unfolding transition in the presence of external mechanical force is likely to be 
first order in all dimensions [TH]. Indeed, molecular simulations of forced-unfolding of proteins 
exhibits a first order phase transition in three dimensions [TUl HI] . However in two dimensions, 
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extensive Monte Carlo simulations performed on a self avoiding walk model in poor solvent [2"U] 
suggests that the folded to unfolded transition could have the hall marks of a second order, in 
accord with the present calculations. Furthermore, scaling analysis shows a second order phase 
transition at a critical force f c [2TJ [22]. Renormalization group analysis of polymer unfolding 
in special lattices also reveals a change in the nature of the phase transition as the spatial 
dimension exceeds two [23J. Thus, the upper critical dimension for force-induced unfolding of 
self-attracting polymers and hence proteins is likely to be three. 

Interestingly, for a range of / < f c {T) the ZL model exhibits a reentrant phase transition 
in which the disordered state is stable at low and high temperatures. As a consequence, f c {T) 
increases as T increases in the low T regime. Force-induced reentrant behavior may be a common 
feature in other model polymeric systems [2TJ[2U[25]. A recent simulation study in the constant 
force ensemble [25J showed that the temperature-dependent force-induced desorption of two- 
dimensional self- avoiding polymer is very similar in shape to the curve plotted in Fig. 2. The 
generality of reentrant phase transitions induced by force, which is similar to cold denaturation 
in proteins, remains to established. It would be of particular interest to characterize such 
transitions, if they exist, in proteins. This would require performing pulling experiments as a 
function of temperature and force. 

In addition to adding to the collection of exactly solvable polymer models exhibiting phase 
transitions, the present work could be used to test approximate theories used to obtain fold- 
ing thermodynamics of proteins using single molecule pulling experiments, which assume that 
extension is an excellent parameter. 
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Figure Captions 



Fig. 1: Conformation of a polymer chain in the ZL model. Two contiguous vertical lines 
define a fold. There is a energy gain per each fold that is achieved by paying a bending penalty. 
The length of each fold is q. A force / is applied to one end of the chain. The calculations are 
done by assuming that the force is along the x axis. 

Fig. 2: Phase diagram in the f — T plane with q = 2 and u/ J — 0.01. At low temperatures 
the critical force increases as T increases reflecting the observed reentrant behavior. 
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